System for synthesizing seismograms from well logs



March 15, 1966 R. A. PETERSON 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS 18am 14am Jamo Mina AyMa/VD A. PErEQsvN,

IN VEN TOR.

r rae/VE y.

March 15, 1966 R, A, PETERSON 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS Original Filed July 5, 1955 14 Sheets-Sheet 2 /Jy/wo/vn A PE TEESU/ IN VEN TOR.

f BY @gf/2,41

March 15, 1966 R. A. PETERSON 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS Original Filed July 5. 1955 14 Shee'cs-Sheei'l 3 Ay/wcwn A. PE Teresa/v,

IN V EN TOR.

By fw March 15, 1966 R, A, PETERSON 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS N VEN TOR. y j

Arran/IEM March 15, 1966 R. A. PETERSON SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS Original Filed July 5, 1955 14 Sheets-Sheet 5 March 15, 1966 R, A, PETERSON 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS Original Filed July 5, 1955 14 Sheets-Sheet 6 .1f/zeig.

www@ w AyMa/vo A .PE rE/eso/v,

1N VEN TOR.

BY M

rraQ/VEM March 15, 1966 R. A. PETERSON SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS 14 Sheets-Sheet '7 Original Filed July 5, 1955 AyMoA/D A. Pamesa/V,

INVENTOR i rraelJEy March 15, 1966 R. A. Pr-:Tr-:RsoN

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS 14 Sheets-Sheet 9 Original Filed July 5, 1955 iff@ f1.

FAZ/0 Jir March 15, 1966 R. A. PETERSON 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS Original Filed July 5, 1955 14 Sheets-Sheet 9 fr@ .f5

JAyMQ/VD A. PETE/@som I N VEN TUR.

BY j/r rrapwey March E5, 19%@ R. A. PETERsoN 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS Original Filed July 5, 1955 14 Sheets-Sheet 10 0 jira /7 DI 6 72 1 70 2:. 2@ e if 'I HQE-l nl.; o)

y l l l J iff@ J9.

.Igny/10H0 Peresm IN VEN TOR.

March 15, 1966 R. A. PETERsoN 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS 14 Sheets-Sheet 11 Original Filed July 5, 1955 IN V EN TOR.

TfaQA/El.

I6: Pa.

R M DA. Peremsa f 336, Ay 0 M March i5, 1966 R. A. PETERSON 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOG'RMS FROM WELL LOGS 14 Sheets-Sheet 12 Original Filed July 5, 1955 lady/wm@ PETEESOM IN VEN TOR.

rra/Q/vey.

March 15, 1966 R, A, PETERSQN 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS Original Filed July 5, 1955 14 Sheets-Sheet 15 jira 3l,

AyMa/vp A, PETEesa/V,

1N VEN TOR. Y @wf/w ,47m/@veg March 15, 1966 R. A. PETERSON 3,241,102

SYSTEM FOR SYNTHESIZING SEISMOGRAMS FROM WELL LOGS Original Filed July 5, 1955 14 Sheets-Sheet 14.

United States Patent Otilice 3,241,102 Patented Mar. 15, 1966 3,241,141@ SYSTEM FR SYNTHESlZiNG SEHSMGGRAMS FRM WELL LOGS Raymond A. Peterson, Altadena, Calif., assigner, by mesne assignments, to United Geophysicai Corporation, Pasadena, Caiif., a corporation of Delaware Original application .filly 5, 1955, Ser. No. 519,801, now Patent No. 3,011,582, dated Dec. 5, 1961. Divided and this application Apr. 3, 1961, Ser. No. 125,020

24 Claims. (Ci. 349-18) This application is a division of copending patent application Ser. No. 519,801, filed July 5, 1955, now Patent No. 3,011,582.

This application relates to improvements in geophysical prospecting systems and particularly to improved systems for use in connection with the interpretation of well logs and more particularly to improved systems for correlating data obtained in various kinds of Well logging systems with the data obtained by means of reflection seismograph systems.

In Ieection sesmograph surveying, as commonly practiced, seismic waves are generated at various seismic Wave generating stations adjacent the surface of the earth. Such waves may be generated at each of such positions, for example, by detonation of a charge of explosive at a shot point located at the bottom of a shallow borehole. The waves so generated at each station travel downwardly through the earth. Upon encountering various subterranean strata or other discontinuities in the formations beneath the surface of the earth the waves are partially reflected and are returned to the surface. There the waves are detected by means of seismic Wave detectors located at various seismic Wave receiver stations and records are made of the returned waves. These records are generally in the form of multiple trace oscillograms in which each trace represents the motion of the earth at a corresponding receiving point, the displacement of the trace being displayed as a function of the seismic wave travel. The depth and the strike and the dip of strata that have reflected the waves may be determined by well known methods from (l) the times required for the waves to travel from each of the seismic wave generator stations to the strata from which the waves are reflected and then to the seismic Wave receiver stations, (2) the geometrical arrangement and relative spacing of the seismic wave receiver stations and the shot point, and (3) the velocity with which the seismic Waves travel. The data so obtained are employed to ascertain the nature and structure of the subterranean formations. Knowledge of such formations is employed to determine Where deposits of petroleum are most likely to be found in the area surveyed. Ordinarily, in reflection seismic surveying, when a survey is made even though the contours representing relative depths of various parts of subsurface formations may be accurately known, it is not at all easy to correlate individual waves recorded on a seismogram with specific subsurface formations. Furthermore, when the character of an individual reflection recorded on a seismogram varies throughout the area under investigation, it is very difficult to ascertain the nature of the corresponding change that has occurred in the characteristics of the formations that have produced the reflection. It is, therefore, very desirable to produce improved methods of correlating the individual formations and changes in the character of such formations with individual seismic waves and with changes in the character of such seismic waves.

This invention relates particularly to systems for correlating logs of wells with seismograms produced in areas in the neighborhoods of the wells. As explained hereinafter, best results are obtained by employing the log of the logarithm ofacoustic impedance as a function of seismic wave travel-time. However, the invention can also be practiced by employing interval-velocity or even the resistivity and even other formation characteristics as the variable in the log.

According to this invention, a system is provided for converting a well log into a seismic wave trace. More specifically in accordance with this invention, after certain preliminary treatment or conversion, a well log is filtered by feeding it through a pulse generator to produce a seismogram trace of the type that would be produced in making a record of seismic waves with a reflection seismograph at the location of the well from which the log was obtained. In the best mode of practicing the invention, well log data are incorporated in a phonographically reproducible well log such as a well log of the variable area type that is a log in which the width of the log varies along the length of the log in accordance with variations in characteristics of the earth in the neighborhood of the well along the length of the Well. For convenience, the record so produced is called a synthesized or synthetic seismogram. By comparing synthetic seismograms produced by this system with actual field seismograms obtained in the same area, it is :possible to identify specific seismic waves recorded in eld seismograms as reliections from specific readily identifiable strata or groups of strata.

Though the invention is described hereinafter particularly with reference to the production of synthetic seismograms from acoustic impedance logs, it has been found that even when other logs are subjected to treatment in accordance with this invention they too may be converted into synthetic seismograms which bear a marked resemblance to actual seismograms obtained in the same area. This similarity probably arises from the fact that there is often a close correlation between the acoustic impedance of various formations with other characteristics of such formations.

In accordance with this invention, interval-velocity logs and other logs may be altered systematically in the laboratory, such as by removal or insertion of sections therein in order to ascertain what effect the truncation or intrusion of various strata may have on actual seismograms obtained in an area under investigation. Furthermore, by way of example, interval-velocity logs may be altered systematically such as by arbitrarily changing velocities at specific depths on the logs in order to determine what effect lithologic changes in a stratum or in a series of strata may have on seismograms produced at various points throughout an area. In these ways stratigraphic changes may be ascertained without the necessity of drilling large numbers of wells over the area under investigation.

The novel features that are considered characteristic of the invention are set forth in the appended claims. The invention itself and the principles upon which it is based, however, both as to its organization and its method of operation, together with various features thereof, will best be understood from the following operation of the principles of the invention and the following description of a specific embodiment thereof when read in connection with the accompanying drawings in which:

FIGURE 1 is a graph Showing how seismic wave travel-time and interval velocity over large depth intervals vary with depth;

FIG. 2 represents Well logs of interval-velocity and resistivity;

FIG. 3 shows a schematic embodiment of the invention;

FIG. 4 is a diagram illustrating schematically how a field seismogram would be produced by seismic wave reections at discreet interfaces;

FIG. 5 is a diagram schematically illustrating how a eld seismogram would be produced Iby reflections from a continuously varying underground formation;

FIG. 6 is a graph illustrating the degree of similarity between the logarithmic approximation and the actual reilection coecient;

FIG. 7 is a diagram employed to explain the principles underlying the invention;

FIG. 8 is a schematic diagram of an embodiment of the invention;

FIGS. 8a and 8b represent logarithmic slits;

FIG. 9 illustrates a synthetic Seismogram produced by the apparatus of FIG. 8, together with a well log from which it was produced and a eld seismogram made in the area of the well;

FIG. 10 is a diagram showing a system for recording a primary seismic wave pulse and a field seismogram;

FIG. 1l is a graph of the relative frequency response :haracteristics of the seismic pulse network and the remainder of the seismograph system;

FIG. 12a is a series of curves illustrating typical outnut waveforms obtained from seismic pulse networks;

FIG. 12b is a pair of curves illustrating two theoretical :ypes of seismic waveforms occurring in viscous media;

FIG. 13 is a block diagram of a seismic wave pulse ietwork;

FIG. 14 schematically illustrates a system for autonatically Iproducing film strip graphs during a logging peration;

FIG. 15 illustrates schematically a system for synthe- `izing seismograms from both interval-velocity and density lata;

FIG. 16 illustrates schematically another system for ynthesizing seismograms `from both interval-velocity and lensity data;

FIG. 17 is a schematic diagram of a system for con- `erting a well log of a characteristic into one repreenting a log of the derivative of the characteristic;

FIG. 18 is a schematic diagram of part of a system for ynthesizing a seismogram from a derivative log;

FIG. 19 is a schematic diagram of part of a system .tilizing a waveform aperture;

FIG. 20 illustrates a variable density waveform slit /hich may -be employed in practicing the invention;

FIG. 21 illustrates a variable area waveform slit which tay be employed in practicing the invention;

FIG. 22 is a schematic diagram of a magnetic tape pparatus suitable for use in the invention;

FIGS. 23a and 23b are diagrams used in explaining ow the magnetic playback heads of FIG. 22 are adjusted;

FIG. 24 is a schematic diagram of a multiple slit ystem;

FIG. 25 illustrates how parts of seismograms synthezed from electrical resistivity and interval-velocity data orrelate with an actual field seismogram;

FIG. 26 is a schematic diagram of a mechanical corzlator;

FIG. 27 is a diagram of a series of logs employed to xplain certain features of the invention;

FIG. 28 is a diagram illustrating the close correlation t a synthetic seismogram with a well log;

FIGS. 29 and 30 are schematic diagrams illustrating 1e method of employing this invention in making vercal velocity surveys; and

FIG. 31 is a graph of interval time versus depth.

General description of seismogram synthesis system This invention is based upon the discovery that certain lnctional operations can be performed on data obtained om well logging operations and the resulting trans- )rmed information may be used as the basis for syniesizing reflection seismograms.

In one method employed to practice the invention, an iterval-velocity log is rst produced in which interval :locities at various depths are displayed as a linear function of total seismic wave `travel-time. It is to be noted that such a log would represent a log of acoustic impedance if the density of the earth was uniform throughout the depth of the Well.

In order to produce such an interval-velocity versus travel-time log, two preliminary logs are first obtained. One is a log of travel-time versus depth, the other is a log of interval-velocity versus depth. From these two logs, or if preferred the log data from which they are produced, the desired interval-velocity versus travel-time log is produced.

Preliminary logs for use in producing an interval-velocity versus travel-time log are illustrated in FIGURE 1. Referring to this figure, it will be noted that the log f L1 represents a plot of travel-time Versus depth, traveltime being plotted as absci'ssae and depths being plotted as ordinates. Likewise, log L2 represents a plot of interval-velocity versus depth, interval-velocities lbeing plotted as abscissae and depths being plotted as ordinates. The same ordinate scales are employed in lboth logs, so that it is a simple matter to correlate the interval velocity at any depth with the time required for seismic waves to travel to that depth.

The travel-time versus depth log L1 may be obtained in numerous ways. For example, a series of char-ges of explosives may -be detonated in a shot hole adjacent the top of a well and the time required for waves to travel through the surrounding formations to various depths in the well may Ibe determined by locating a hydrophone at the various depths in the well as each shot is taken and measuring the time required for the Waves to travel from the respective shots to the hydrophone while at each of the corresponding depths.

Likewise, continuous interval-velocity logs may be determined in numerous ways. In one method, for eX- ample, seismic Waves may be successively generated at various depths in a well and the difference in the time required for the waves to travel to hydrophones maintained at xed distances from the source of the waves may be measured at the different depths. The interval velocity is computed by dividing the inter-hydrophone distance by the suitably corrected interval time at each depth. Such a system has been disclosed in detail and claimed in copending United States patent application Ser. No. 462,062, filed Oct. 13, 1954, now Patent No. 3,170,136. Other systems for determining interval-velocity logs have Ibeen described by Summers and Broding in an article entitled Continuous Velocity Logging, published in Geophysics, vol. XVII, No. 3, July 1952, pages 598 to 614, and another by C. B. Vogel in an article entitled A Seismic Velocity Logging Method, published in the same issue at pages 586 to 597.

In either case, the travel-times so determined and the continuous interval velocities so determined may be tabulated as a function of depth or they may be plotted as a function of depth as illustrated in FIG. l.

In practice, interval-velocity measurements should be made at different depths over small intervals of distance less than about 20 feet or l0 feet in order to be most useful in this invention. Thus, wheras the interval-vel ocity log L2 illustrated in FIG. 1 shows measurements over intervals of at least feet, a much more detailed log is normally required for use in this invention. In FIG. 2, there is illustrated an interval-velocity log L3, in which interval velocity over small intervals of about 10 feet is shown plotted as a function of depth. In the same Iigure, a resistivity log L4 for the same well is plotted over the lsame depth range. It will be noted that there is a striking similarity between these logs. Advantage of this fact is taken in accordance with this invention to produce synthetic seismograms from resistivity logs as well as from interval-velocity logs.

In preparing an interval-velocity log to use for producing a synthetic seismogram, a phonographically reproducible Well log in the form of a variable area film strip log L5, as illustrated in FIG.3, is prepared. As will be more fully explained hereinafter, in the best logs useful in this invention, the lateral width variable area, or Opaque part, of the record L is -rnade proportional to the logarithm of the interval velocity rather than proportional to the interval velocity itself, while the coordinate along the length of the variable area log is made proportional to the total corresponding seismic wave traveltime rather than proportional to the depth. Thus, in the variable area log of FIG. 5, the width w of the opaque part of the variable area record at any point along its length is proportional to the logarithm of the interval velocity at the depth to which seismic waves travel from the shot point to a receiver in a time T. Such a log is readily prepared from the data of the log L1 of FIG. 1 and the data of the log L3 of FIG. 2. A method for producing such a log automatically while a well is being logged is described in detail hereinafter.

In accordance with this invention, the variable area log L5 is fed into a specially designed seismic wave synthesizer or transducer M, where the variable area log L5 is scanned to produce a series of electrical signals which are combined to produce a synthetic seismogram S1 as illustrated generally in FIG. 3.

As more fully explained hereinafter, the variable area film strip log L5 is scanned in a phonographic reproducer operated at constant speed. The output from the reproducer is filtered by a filter which simulates the action of a seismic wavelet propagated into the subsurface formations by the explosion of the dynamite charge during the making of a reflection seismogram in the field. The resulting output signal from the electrical iilter network is recorded as a function of time and becomes the synthesized or predicted seismogram S1 shown in FIG. 3 in which the trace represents the manner in which the velocity of the ground varies as a function of travel-time at the point where the reflected waves are received. In FIG. 3, the interval-velocity log L5 and the synthetic seismogram S1 are both shown as functions of wave travel time and are 1placed side-by-side wit-h time increasing from the top of both records to the bottom thereof to facilitate visual correlation of corresponding parts of the records. It will be understood that the oneway travel-time of the waves from the shot point to the various strata is one-half of the two-way travel-time from the shot point to a reflecting stratum and return to the shot point. This relationship is indicated in FIG. 3, where the time interval t on the syntheic seismogram S1 representing two-way reflection time occupies the same length of record as the time interval labeled 1/zt corresponding to one-way travel-time on the well log L5.

For convenience hereinafter a well log exhibiting a well formation characteristic plotted as a function of travel time is sometimes referred to as a travel-time well log, and one which exhibits the well formation characteristic plotted as a function of depth is sometimes referred to as a depth well log.

It is significant in the method of this invention that a synthesized or predicted reflection seismogram is produced which is free from all surface disturbances, multiple refiections, extraneous noises, etc. The seismogram is, in effect, a true picture of a simple reflection process which is to be expected from the nature of the physical characteristics of the formations through which the borehole was drilled.

General theory of seismogram synthesis It is well known that a plane compressional elastic wave incident normally upon an interface exhibiting changes of seismic-wave velocity characteristics or density, or both, will be partially reflected therefrom. In fact, under these ideal conditions it is possible to cal- Gg culate the reflection coeillcient which, when multiplied by the original incident wave function, determines the amplitude of the reflected wave relative to the incident wave. The remaining portion of the original wave passes through the interface and is propagated further. More particularly, if the incident Wave is originally traveling in a medium having an acoustic impedance Z1 defined as the product of the density p1 and the acoustic velocity of propagation c1, and if the other medium just on the opposite side of the interface is defined by corresponding symbols having the subscript 2, then the reflection coefficient RX is defined by It should be observed that the algebraic sign of the reflection coefficient RX is determined by the relative magnitudes of the acoustic impedances Z1 and Z2 adjacent the interface.

A diagrammatic representation of the physical process of the reflection of seismic waves at acoustic interfaces is shown in FIG. 4. A combination density-velocity, or acoustic impedance, log exhibiting distinct contrasts at three different depths is shown at the left of the figure. The horizontal coordinate Z of this log is proportional to the product of the density p and the characteristic velocity c of the formation. This combination of the elastic properties of the formation is plotted as a linear function of vertical travel-time as shown. It is assumed that a simple shot pulse in the form-of a primary seismic pulse or wavelet 1 enters the medium at point 2 of the surface 6. As the incident Wavelet 1 propagates downwardly, it strikes the interface 3 of the densityvelocity characteristic at time 1/23, thus being partially reflected and partially transmitted. The reflected wavelet 3 exhibits the same polarity but has a different amplitude than the incident wavelet 1. The remaining portion of the incident wavelet 1 continues to travel past `interface 3 and encounters interface 4 at time 1/2t4. Again a fraction of the incident wavelet is reflected upwardly, thus forming reflected wavelet 4. Similarly, a further fraction of the remaining wavelet is reflected upwardly at interface 5 at time 1/2t5, thus forming reflected Wavelet 5. In this case, however, interface 5 exhibits a decrease in acoustic impedance and therefore the reflected wavelet 5 is inverted.

To schematically illustrate the resultant seismogram, a horizontal time scale is indicated commencing at the time to of the initiation of the incident wavelet 1. Each of the reflected wavelets is plotted at a position on the scale corresponding with the time of arrival of the wavelet at the surface. In this way, a simple seismogram 7 is formed showing the incident wavelet 1 occurring at time to and each of the reflected wavelets 3', 4, and 5 arriving at the surface at timesv t3, t4, and t5 respectively.

Several points may be observed from this simple illustration. In accordance with Equation 1, the amplitudes of the reflected waves are proportional to a function of the change in acoustic impedance Z of the density-velocity log. Secondly, the polarity of the reflected Wave is determined by the direction of change of the acoustic impedance Z. Note that the interfaces 3 and 4 represent an increase in the value of Z, While interface 5 represents a decrease. Consequently, reflections 3' and 4' each begin with an initial downward swing while reflection S begins with an initial upward swing. And thirdly, the time required for a wave to travel from the shot point to a reflecting interface and return is twice the traveltime required for a wave to travel from the shot point to the reflecting interface. The former travel-time is measured in a reflection seismic prospecting system while the latter is measured in a well shooting system. It is to be noted that the interfaces shown in FIG. 4 may be designated generally by the integer j, and that the formation directly above the interface j has an acoustic impedance Z5, while the formation directly beneath the interface j has in a-coustic impedance Zj+1, and that the wavelet reflected rom the jth interface is j.

The density-velocity log illustrated in FIG. 4 is, of zourse, highly simplified. In practice, the distribution of he values of the acoustic impedance with depth is typi- :ally quite irregular, as schematically illustrated in FIG. In view of the complexity of this illustration, one nay ask the general question, At what time are seismic eflections to be expected on seismograms recorded at he well location, as predicted from the observed variaions of velocity and/or density with depth at that locaion? Simple inspection of the density-velocity log re- 'eals that the answer to this question is quite complicated necause of the highly irregular distribution of acoustic mpedance. However, by means of this invention, prelicted or synthesized seismograms can be produced diectly from the density-velocity, or acoustic impedance, ogs.

In accordance with the present invention, the density- 'elocity log of FIG. 5 is considered to be made up of a arge number of small discrete discontinuities or steps 10. lach of the steps 10 exhibits a definite reflection coeffiient and gives rise to a small reflection wavelet which eturns to the surface 6. The total reflected wave on the esultant seismogram 7 is then made up of the linear ummation of all of the reflected wavelets.

Further, in accordance with the invention, the density- 'elocity versus depth log of a well is transformed into a few form in which the logarithm of the density-velocity 'ariations is plotted as a function of the total seismic vave travel-time instead of the depth. This transformed 3g is then converted to a proportional electrical signal in /hich the amplitude of a voltage or current varies as a unction of time in the same way that the logarithm of 11e impedance varies as a function of travel-time. It is nderstood, of course, that the time scales of the voltage vave need not be the same as the travel-time scale. Howver, unless stated otherwise, these two scales are assumed 3 be the same in the description that follows. The sigal produced is applied to the input of a seismic-pulse r synthesizing network in the synthesizer M. The haracteristics of the network are adjusted so that the :ansient produced at the output of that network when a tep function signal is applied to the input has the form f the primary seismic pulse which is traveling downvardly through the formations. In this way, the stepvise characteristic of the transformed density-velocity log cts as a series of successive input step functions of varyng size, polarity, and spacing. The output of the net- Jork is then made up of the linear summation of all of 1e reflection wavelets applied at the proper times and lith the proper phases and amplitudes. Thus the desired ynthesized seismogram is obtained.

This process has been diagrammed schematically in "IG. 3, which shows the transformed acoustic well log as a function of time) as applied as an input signal to he seismic pulse network. The resulting output is the esired record, or synthetic seismogram, representing the um of reflected seismic waves arriving at the surface. n FIG. 3, the output signal has been shifted by the mount of time necessary to correct for the time delay in- :oduced by the seismic pulse network. As will be seen iter, by shifting a field seismogram relative to a travelime well log by such a time delay, the correlation of the eld seismogram with the well log is facilitated. In this Jay the close correlation between the general waveforms f the input and output signals is brought to light.

Mathematical analysis A more detailed explanation of the mathematical theory lnderlying the invention will now be made since this .nalysis will aid in understanding the invention. In carryng out the analysis certain simplifying assumptions are nade, even though the general results of the theory apply o under a wide variety of conditions. Conditions assumed to exist solely for the purpose of developing a simple theory are:

(l) The media are considered to be horizontally stratified and to exhibit no changes in elastic properties in the horizontal direction.

(2) Multiple reflection phenomena are disregarded.

(3) The amplitude and shape of the incident seismic wavelet are considered to be invariant with depth.

(4) The incident wave is considered to be plane and infinite in extent, i.e., spherical divergance of the wave is neglected.

In practice, assumptions (l) and (4) are closely approximated in many areas. The effects of attenuation of seismic waves with depth (assumption (3)) are effectively overcome in practice through the use of seismograph amplifiers employing automatic gain control systems such as automatic Volume control systems and expander" gain controlling systems and combinations thereof. Account of changes in shape may be taken in practice by Varying the characteristics of the seismic pulse network as a function of travel-time. Absence of the effects of multiple reflection phenomena will be shown to be very advantageous.

By employing the expression for the reflection coeflcient RX given by Equation l, the process of seismogram formation illustrated schematically in FIGS. 4 and 5 may be defined in mathematical terms. Assume, with respect to FIG. 4, that an arbitrary plane pressure wave Pi(t) is propagated vertically downward from the surface thus striking the interfaces at normal incidence. The incident wave P1(t) travels at a velocity c and arrives at the first acoustic impedance interface 3 after an appropriate travel-time 1/2t3. Some of the incident energy travels through the interface 3 while a fraction is reflected upwardly by the interface. More particularly, in accordance with Equation l the reflected pressure wave arriving at the surface is given by In this equation Z3 and Z4 are the values of the acoustic impedance on the upper and lower sides of the interface 3 respectively as shown in FIG. 4. Time t3 is the total travel-time of the initial part of the wave from the surface 6 to the interface 3 and back again to the surface. It is seen that the quantity within the brackets is the reflection coefl'icient RX for interface 3. The reflected wave represented by Equation 2 is propagated upwardly, ultimately arriving at the surface 6 where it is recorded as wavelet 3 on the seismogram 7.

The remaining part of the incident energy which passed through interface 3 continues to travel downwardly and ultimately encounters the interface 4 where a further fraction of the incident energy is reflected upwardly, the remainder traveling through interface 4 in a downward direction. The reflected wavelet from interface 4 has the same form as that reflected from interface 3 and is given by Equation 2 where t3 would assume a new value t4 and the subscripts 3 and 4 of the acoustic impedance terms Z would be replaced respectively by the subscripts 4 and 5. Likewise, the `energy which passes through interface 4 continues to travel downwardly ultimately to encounter interface 5 where the reflection process is again repeated. At the latter interface it should be noted that the acoustic impedance below the interface has a lower value than that above the interface and hence from Equations 1 or 2, the sign of the reflection coeficient RX becomes negative. This negative sign indicates that the wavelet 5 reflected from interface 5 is inverted in phase or polarity as compared to the incident wave, all as shown schematically in FIG. 4.

An expression of the form of Equation 2 may be written with additional terms and suitable time delays between the terms to express the sum total of the reflected Wave information received at the surface 6. Such an expression will have the appearance of Equation 3,

where the various symbols have the meanings indicated in Equation 2. It will be understood that the expression P10) represents the manner in which the pressure of the down-going wave varies as a function of time and that, generally speaking, the expression P(z-tj) represents the manner in which the pressure of the down-going wave varies at the jth interface. All of the expressions P(t), Pi(t-t3), P10-t4), etc., have the same wave form, but occur at different times. The down-going wave incident on the jth interface is P10-tj). This wave has the same form as the original way Pi(t), but is delayed by a time tj, the time required for the original wave to travel from point Z to the corresponding jth interface, and from there to the surface. In Equation 3, the expression PTU) represents the manner in which the pressure of the upwardly-traveling train of reflected waves that is received at a point near the surface of the earth varies as a function of time. In all of these cases, the expressions PU), Pi-tj), and Pr(), represent the waves arriving at the surface, the Itime t being the total time elapsed since the wave traveled downwardly past the point 2. In these expressions, no account is taken of Changes in amplitude of the various waves caused either by reflection, or by other attenuation. It is to be noted that in Equations 2 and 3 when t-tj= (5) where j=3, 4, 5, etc (6) As previously mentioned, the actual variation of acoustic impedance Z as a function of depth is found to be an extremely variable characteristic as illustrated in FIG. 5. It can be seen, therefore, that to carry out the process indicated by Equation 3 for any practical case would be an extremely laborious process indeed.

However, at this point a great simplification is introduced, the discovery of which made the present invention possible. More particularly, in Equation l it is noted that for small differences AZ Abetween the values Z1 and Z2, the reflection coeliicient at the depth x is approximately given by the expression AZ dZ z zz (7) The latter function of Equation 7 is recognized as onehalf the differential of the natural logarithm of the acoustic impedance Z. Thus,

P. e) =d-U] Z)P. c-fa 9) This equation expresses the fact that the wave reflected from an acoustic interface is proportional to one-half of the corresponding incremental change in the value of the logarithm of the acoustic impedance Z.

In using this expression, the earth is visualized as being made up of a large number of rock layers 1I (FIG. 5), very thin compared with the shortest seismic wave lengths of interest, and each having an acoustic impedance that is generally somewhat different from those of the adjacent layers above and below it. The downwardly traveling incident wave or seismic pulse is partially reflected from each interface between successive layers and the seismic wave train received at the surface is the sum of al1 of these small reflections. It should be noted that the logarithmic approximation of Equation 8 affords a significant simplication since it vitiates the need for any step-by step calculation procedures which would otherwise be necessary in determining the values of the reiiection coetiicient corresponding to the multitude of rock layer interfaces.

It is evident that the approximation inherent in the transformation from Equation 2 to Equation 7 becomes better as the size of the discontinuity decreases. In the limit for which the acoustic impedance Z becomes continuous with no step discontinuities whatsoever, the approximation is exact. In dealing with practical well velocity data, however, the curve of the acoustic impedance Z often assumes a discontinuous step-wise characteristic and for this case it is Valuable to know the extent of the approximation involved in Equations 7, 8, and 9. In in order to make this evaluation, consider the following:

The equation is first transformed into a dimensionless form by defining d "n+1 (u) From Equation 8 the approximation for the reflection coefficient RX is written in the form Raz-2 2 2 2 A series expansion for the last term of Equation 8 is:

2 )url-l Inspection of Equation 13 reveals the fact that the first term of the expansion is, in fact, the exact value of the reflection coeliicient RX as shown by Equation 11. All of the higher order terms therefore represent slight errors inherent in the repesentation of the reflection coeiiicient RX by the term ll'l ,u. 2

It is now apparent that the real approximation inherent in the transformation from Equation 2 to Equation 3 is the assumption that the value of lnlr 2 is given by only the rst term of its series expansion.

In order to evaluate the practical limits of this approximation, a curve relating the actual reflection coefficient and approximate reection coetiicient (abseissa) l. i are inherent in the present invention, it is reasonable to extend the range of operation further than that just described. In fact, at a value of reflection coecient of +05 the error is only 10%. Thus it is seen that the assumption that the quantity is valid for a large range of values. The curve of FIG. 6 may be employed to determine the actual range of the equivalents for any desired degree of accuracy.

Having thus shown that Equation 9 accurately represents the characteristics of wavelet reflected from any particular interface of acoustic impedance, it is now desired to continue further with a simple mathematical analysis in order to bring out clearly the principles underlying the invention.

FIG. 7 illustrates a transformed strip graph for which the abscissa coordinate is the total seismic travel-time t, while the ordinate is proportional to the logarithm of the acoustic impedance Z. This form of graph can be obtained in different ways, but fundamentally the depth scale of the well log is converted into a total vertical travel-time scale by adding the small time intervals At corresponding with each of the interval velocity values measured in the well and multiplying the result by two in order to account for travel in both directions.

In FIG. 7 the curve ln Z(t) is shown approximated by a series of step functions. The rst is applied at time t=Ar. The second is applied at time t=2Ar and is of a height necessary to reach the curve ln Z(t) at that point. Thereafter, appropriately sized step functions are added at equal intervals of time AT in order to approximate the curve. By making the size of A-r sufficiently small, the approximation to the curve ln Z(t) produced by the summation of step functions may be made arbitrarily close to the actual curve. In effect, each of the small step functions approximating the curve acts as an acoustic interface exhibiting a reflection coefficient proportional to the size of the step. Thus, each f these small steps produces a small reflection when the incident pressure wave Pi(t) is propagated past it. By plotting the curve of FIG. 7 in terms of a total travel-time instead of in terms of depth, the ensuing analysis becomes highly simplied.

More particularly, in response to an initial applied plane pressure wave having a form Pi(t), the curve of FIG. 7 produces a reflected wave train having the form given by A-r Af +ATPitt-2AT) [AUALTZ +ATP, (t-naf) [MLTZLAT (14) where A (l-AI-l-Z-l] =the ratio of the incremental change in T nAr ln Z to AT over the nth Ar interval (l) This quantity approaches the value of the derivative QLLIH'Zw (16) 1s Af becomes vanishingly small. Now let the quantity HATIT (17) l`hen if there is a limit it is given by P,(t) :lAoPiG-T) ln Z(T)A1 (18) Cil However, since ln Z(t) actually varies continuously and also since the value of A1- may be made arbitrarily small, the final expression becomes If any discontinuities occur at any point these may be taken into account by well known mathematical methods. This equation expresses the fact that when an initial pressure wave Pi(t) is propagated downwardly through the changes of acoustic impedance Z expressed by FIG. 7, the resulting reflected wave is an integral representing the continuous runningr cross-correlation between the derivative of the logarithm of the acoustic impedance Z as a function of time and the initial applied pressure wave P(t). It is to b-e noted that Equation 19 specifies only the reflected wave and does not include the seismic wave pulse that would travel directly to a geophone located very near the shot point. Equation 19 is in the form of a linear superposition integral. Moreover, the Expression 19 represents the voltage that would be -observed at the output terminals of an electrical network whose step function response is Pi(z) and which is driven by an input waveform having the characteristics shown in FIG.

7. Such a network automatically produces the integrak tion indicated in Equation 19 and hence produces a record of seismic waves that would be reflected by a series of strata in which the acoustic impedance is varied in accordance with the function Z. Use is made of the discovery of `this relationship in the synthesizing of seismic waves.

In one form of this invention, an electrical waveform which is proportional to the logarithm of the acoustic impedance as a function of the total seismic travel time t is formed, and this waveform is passed through a network whose response to a step function has the same waveform Pi(t) as the incident seismic pulse. The output from the network represents the instantaneous particle velocity of ground motion and is recorded by an oscillograph camera to produce the desired seismogram trace. In addition, this output wave is applied to a circuit that is the electrical equivalent of a geophone and a seismic amplier and the resulting signal also is recorded in the form of a seismogram trace.

Electrical filtering system In FIG. 8 there is illustrated a practical `system which is particularly adapted to carrying out the process of synthesizing reflection seismograms in accordance with this invention. As previously mentioned, it is desirable to generate a signal which is instaneously proportional to the logarithm of the acoustic impedance of the formation surrounding the borehole as a function of the total vertical seismic wave travel-time. This process is accomplished in the system of FIG. 8 by means of the photoelectric scanner 50. This device scans or reads the instantaneous amplitude of the opaque portion of the variable area film strip graph 52 as the latter is rotated at a constant speed past a light source and scanning slit. The film strip graph 52 exhibits the characteristic of the curve of FIG. 7, the relative height transverse the length of the opaque portion at any instant being proportional to the logarithm of the acoustic impedance Z of each particular stratum, and the coordinate along the length or periphery of `the film 52 being proportional to the total travel-time of the seismic waves traversing the path from a shot point at the surface of the earth down to the particular strata and back to the surface again. Also shown on the strip graph 52 are step functions 54 which have accurately measured heights and are employed for calibration purposes as hereinafter described.

The photoelectric scanner 5t) is constructed along the general lines of a device described by Washburn and Olesen in an article entitled Transient Testing of Seismic 

1. IN A SEISMOGRAM SYNTHESIZER: CROSS-CORRELATION MEANS HAVING AN INPUT AND AN OUTPUT, SAID CROSS-CORRELATION MEANS HAVING SUCH A CHARACTERISTIC THAT A VARIABLE INPUT SIGNAL APPLIED AT A POINT IN SAID INPUT GENERATES AT SAID OUTPUT A COMPOSITE SIGNAL WHICH IS THE CROSS-CORRELATION FUNCTION OF A SESIMIC WAVE PULSE AND A FUNCTION OF SAID INPUT SIGNAL, MEANS FOR CREATING AN INPUT SIGNAL WHICH VARIES WITH 